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CALCULATION OF ENERGY DEPOSITION DISTRIBUTIONS 
FOR SIMPLE GEOMETRIES 

I. INTRODUCTION 


Parallel plate ionization chambers are to be used ( 1 ] as part of an instrument to 
measure the charge spectrum of heavy nuclei in the galactic cosmic-ray flux. The ionization 
chamber measurement determines the energy deposited by the cosmic ray along its 
trajectory through the chamber while other portions of the instrument determine the 
distance traveled through the chamber. (The trajectory is approximately a straight line 
unless the cosmic ray has a nuclear interaction.) From these two measurements, the energy 
deposited in the chamber per unit pathlength of the cosmic ray can be determined. The 
stopping power or mean energy lost per unit pathlength in material - dT/dx is directly 
related to the incident particle’s charge. For nuclei up to about iron, calculations of 
stopping power based on the first Born approximation appear to be valid |2,3). These 
calculations predict that the stopping power is proportional to the nuclear charge squared. 
For nuclei above iron, higher terms in the perturbation theory series must be taken or a 
partial wave expansion must be used (3,41 to derive the charge-stopping power relationship. 

However, the stopping power gives us only the a"erage energy lost by the particle in 
the detector. A beam of particles of the same trajectory, charge, and energy will produce a 
distribution of energy losses in the ionization chamber. This fluctuation in energy loss in a 
thin absorber, sometimes called “straggling,” is the result of the finite number of events that 
occur in the ionization energy loss process. Energy transfer to the electrons in the absorber 
material produces excitation and ionization. Many of the electrons gain considerable energy. 
The probability of large energy transfers, which depends upon the impact parameter in the 
classical picture, becomes smaller with higher secondary electron energy so that the 
fluctuations in the number of events become greater at higher energies. As a result, the 
energy loss distribution for fast charged particles in thin absorbers has a considerable width 
and is skewed to the higlr energy side. The problem of calculating this distribution was first 
attacked theoretically by Landau (51 for thin absorbers where only a small fraction of the 
particle’s initial energy is lost. He used Laplace transforms to solve the continuity equation 
of the distribution function of energy loss. In addition to the assumption that the energy 
loss was small, his derivation places no upper limit on the energy transferred to a single 
electron of the absorber. This approximation limited the validity of his solution to 
electrons, positrons, and heavy particles of very high energy or in very thin absorbers. 
Symon (61 solved the problem more generally by using numerical techniques both for thin 
absorbers where the average energy loss was less th.^n 10 percent of the initial energy and for 
thick absorbers where the average energy los-s is more than 10 percent of the initial energy. 
Vavilov (7] followed Landau’s basic approach, but he obtained a more rigorous solution 
without the approximation of unbounded single energy transfers. Thus, his results could be 
applied to a wider range of particles and energies though still restricted to thin absorbers. We 
will use Vavilov's general method of solution extensively in Section 11. 
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For a typical primary cosmic ray with kinetic energy in the GeV/nucleon range, a 
number of energy loss mechanisms are available — scattering of electrons (delta rays), 
excitation of atomic electrons to higher bound-energy states, production of Cerenkov 
radiation (the radiation produced when a particle passes through a dielectric at a velocity 
faster than the phase velocity of light in the medium) {81, and production of 
bremsstrahlung radiation (radiation that is produced when a charged particle undergoes 
acceleration). Production of delta rays and excitation is the most important energy loss 
process contributing to the total stopping power for the particles considered here. Cerenkov 
losses are responsible for only a small fraction of the relativistic rise of the stopping 
power (9| , and bremsstrahlung losses are negligible because of the high mass of the nuclei. 
We will confine ourselves to consideration of the interactions with electrons. 

Cosmic rays can also undergo nuclear interactions with the chamber material. 
Waddington et al. (10) have measured the interaction mean free path for cosmic rays in 
nuclear emulsions. Their results show that the mean free path varies, with minor energy 
dependence (energy range was 100 MeV/nucleon to 30 GeV/nucleon), from about 19 cm of 
emulsion (the equivalent of about 46 g/cm* of air) for alpha particles to about 10 cm of 
emulsion (27 g/cm* of air) for very heavy primaries with nuclear charges greater than 20. 
Because the ionization chambers considered here are relatively thin (x < l.Og/cm*) 
compared with the interaction mean free path, few p''clear interactions will occur, and they 
are identified by the trajectory measurement in the experiment ( 1 1 and treated separately. 

The problem to be attacked here is the result of the subtle distinction between the 
energy lost In a given element of volume along the cosmic ray’s path and the energy 
deposited there, since energy can be transported in and out of the volume element by 
secondary particles and radiation. Thus, energy loss depends on the type of material in the 
volume, whereas energy dcpos-ied depends on the surrounding geometry as well. Since delta 
rays transport practically all energy that is removed from the region where it is lost by the 
primary particle, let us consider how geometry affects their production, transport, and 
energy deposition. For a small region on the path of a heavy cosmic ray through a piece of 
material, the production of a high-energy delta ray will cause a large fluctuation in the 
energy lost in the region; but, since the delta ray will leave the region with little energy loss, 
the fluctuation in energy deposition will not be as large. Conversely, electrons from other 
portions of the path may scatter into the region, contributing to energy deposition in the 
region. More of the energy of lower-energy delta rays is deposited in the region of 
production. Finally, at some deita-ray energy, depending on the size and shape of the 
region, energy loss is equivalent to energy deposition because few of the delta rays escape 
from the region relative to the number produced there. 

F.ven if the relationship between the mean energy deposition and charge was known, 
problems would exist because the production and transport of the delta rays are random 
processes and the energy deposited in a region by identical cosmic rays will fluctuate from 
event to event depending on the number and energies of the delta rays produced. For an 
ionization chamber, this results in an uncertainty in the determination of the exact charge 
of the cosmic ray. Thus, the interesting quantity that we will attempt to calculate here is the 
frequency distribution of energy deposition by a heavy cosmic ray considering the 
fluctuations in delta-ray production and transport. 
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To attack the delta*ray problem, we will use the procedure of dividing it into two 
problems according to the delta-ray energy. The breaking-point energy e, is selected to 
meet two criteria. First, it must be large compared with the binding energy of electrons in 
the material so that electrons scattered at energies above this energy may be treated as free 
electrons for the scattering cross-section calculation. This places a lower bound on ei in 
the 10- to lOO-KeV range (111. Second, since we will assume that all electrons with energies 
below C| deposit their energy in the volume element where they are produced, so that 
energy loss Is equivalent to energy deposition, the distance that electrons with energies 
below 6 1 can travel must be on the order of the smallest dimension of the production 
volume element. For the calculation shown here. ci was 80 or 100 KeV because the 
ionization chambers considered are relatively thin. For other geometries, C| could be much 
higher. 


We will calculate the cnergy-loss/deposition distribution caused by low-energy delta 
rays and excitation using a modification of the Vavilov energy-loss distribution (7) to be 
derived here. The modification consists of introducing a cutoff energy above which all delta 
rays arc ignored. For the delta rays abo " ei , we will calculate the production number and 
energy distribution, and then determine energy deposited in the detector by these electrons 
using Monte-Carlo techniques. Since delta rays produced in the chamber electrodes and 
other parts of the instrument can scatter into the chamber gas, the Monte-Carlo transport 
will consider delta rays produced throughout the instrument. 


II. THE DISTRIBUTION OF ENERGY LOST 
DUE TO LOW ENERGY TRANSFERS 


A. Development of the Restricted Energy>Loss Distribution 

The distribution of energy loss due to low energy transfers is derived in this section. 
Tlie derivation shown here follows Vavilov’s original derivation of the energy-loss 
distribution (71 where all energy transfers were considered. Here, we exclude high-energy 
transfers by the introduction of a cutoff energy C| . As we will show, the introduction of 
the cutoff energy in the derivation is not equivalent to replacing the maximum transferable 
energy in the Vavilov distribution equation by a lower energy. 

If we define f(X,A)dA as the prob.ihility of an energy loss between energy A and 
A + dA because of interactions wh.ch re.suit in losses below ci for particles traveling a 
distance X and w(E,e)dc as the probability per unit length of an energy loss e by the 
particle at energy B , then the continuity equation for f(X.A) is 


8f(X,A) 

ax 


b 

- f a>(E,e)f(X,A - €)d€ - / w(E,c) f (X, A)de 
o o 


(1) 
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[See Appendix E for a derivation of equation (l).l Wc assume that E > A so 
that w(E,e)does not change as a function of E along the particle path, e, is the cutotf 
energy above which the delta rays are ignored. It is assumed that u(E,e) = Ofor e > «i . 
l^us, b is equal to A for A < C| and equal to C| for A > ei . 

To use Laplace transforms on equation (1), we note that 


^(X,P) a P / f(X,A)dA 


^ T <»> 




^(E,P) a p J aj(E,€)de 


C+l« P€,/_ 

= -ti 


(2) 


where C is an arbitrary constant. When b equals A , using the convolution theorem on 
the first term on the right ol equation (1 ), wc obtain 


P) ^ ii-(E.P)<;>(X,P) _ ^(x,P) w(E,e)de (3) 

9X P o 

Solving the differential equation, we get 
«(X,P) = 


0(0, P) exp 


I a>(E,e)dc|x| 


_pg 

J e *• w^,€)d€ 
o 


(4) 
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• 4 ' * ; f-4"' * 


thus. 




0(X, P) = ^(0 , P) exp < 


f w(E, €)de - f w(E, c)de 


,e)de x| 


. (5) 


The incident particles will have no losses; thus, 


f(0,A) = 6(A) 


where 5(A) is the Dirac delta function and 


^(O.P) = P / 6(A)dA 

o 


« P 


By transforming equation (5), we have 


Ol® PA 

“ InT / 

C-|co 


^P exp ^ 


f e”^* w(E, e)de 


, i 


- J w(E,e)de 


XV\dP 


( 6 ) 


I 1 
i 


or 


f(X.A) 


1 

2rri 


C+i« 

r 


C-ico 


exp 



+ X /* w(E,e)(e“^^ - l)de 


dP 


(7) 


For b equal €| , we get the same equation because of the way w(E,c) is defined. This is 
basically Vavilov’s form of the expression except that in his case C| is the maximum 
transferable energy. We may rewrite the exponent in equation (7) as 


I = P 


A - X J eu(£,e)d£ 


ft Pc 

X / w(E,e)(e + 

o 


Pc - l)d€ (8) 


where distant collisions between the particles and electrons play an important role in the 
first term only. Tlic integral in the first term can be written as 

€ € 

max max 

Ij = J cu(E,e)dc - J eu)(E,c)de (9) 

o Cl 


where fmav is the maximum possible energy transferred to an electron given by 


e 

max 


1 + 


mc^ 

2m /niV 
M \M/ 


Here, m is the electron mass, c is the speed of light, 0 is the primary particle velocity 
over c . M is the primary particle mass, and y is defined as 


r 


1 

^/l - /J* 
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(Appendix A shows a derivation of equation (9) is recognized as 

the definition of collision stopping power S(E) . Thus, 


I = P\A - X 


max 


S(E) - / €w(E,e)de 

ei 


+ X /* w(E,e)(e"**^ + Pe - l)dc 


( 10 ) 


If we confine ourselves to heavy energetic charged particles, we may use the results that 
have been obtained for w(E,c) by applying the first-born approximation if ci is assumed 
to be higli enough above the electron binding energy in the first-term integral. The 
probability per unit pathlength for the production of an electron with energy e is given 
by [121 


oj(E,e) 


27tNZV /, 


€ 

max 


( 11 ) 




where N is the number of electrons per unit volume, Z is the atomic number of the 

primary particle, and e is the electron charge. Defining | as s 




2ffNZV 

~mcW~ 


( 12 ) 


and substituting equation (I ! ) into equation ( I) we have 


I = P 


max 


A - XS(E) -if i 1 - 


max 


d€ 


o \ max / 


+ Pe - l)de 


(13) 
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Performing some of the integrals and replacing XS(E) , which is the mean energy loss, 
with A , \'>e obtain 


I = P 


A - A + X5 




max o2 . 


- P* + 

€i ^ e 


max; 


+ X4 


-ClP 


\ e 

* \ VMO 


t / \ 'max/o 


Cl , -eP 
-* 1-e 


/ 


de 


Simplifying slightly and noting that [13] 


€i , -eP 
-* 1-6 


/ dc = t + In ejP - E^(-eiP) 


where f is Euler’s constant 0.577 ... andEj(-ciP) is the exponential integral, we obtain 


I = p|a - A - X5^1n + 1 - 


Letting 


y = €jP 


K = ^ 


8 



A - A 


- kAd + /?*+!-£ 

\ 'max 


and substituting into equation (7). we obtain 


+ ^Y + Y - Ej(-Y)j| jdy 


Choosing C' = 0 for the contour of integration and letting iu = Y , we obtain 


exp K 


</i * 


f(X,A) = 


f exp|l\'t 


- /lu + - E^(-iu)l|i< 

( V ^max/ * * 


Since c'*^ = cos u - i sin u. 


In I u = 


In “^lui = Ijt/ 2 + Inlul u a 0 
In = -lir/2 + Inlul u < 0 







'•.• '•j 


and from Reference 13 


( C (|u|) - i S (u) + i7r/2 u a 0 

Cj(-iu) = / 

I Cj(|u|) - i S^(u) - i7r/2 u < 0 


f(X,A) = 


expKfl + f^j „ 

— J ^l'’ 


K <cos u - I sin u - jiu + j[ln|u| 

* \ max / 

- a(Iu|) + iS^(u)]||du 


where Cj(u) is the cosine integral and Sj(u) is the sine integral. Regrouping into real and 
imaginary parts 


'Ji ■<■ ^ 


f(X,A) = 


u Sj^(u) - cos u > + iXu + 


+ ■ S.(u) + sinutidu 

max ) ; 


j exp I x| [Inlul - a(|ul)] 

-® ( ( max 

u I + i\u + i« |u[ In I u I - C.( | u 
sin u| |du 


Letting 


ei(u) = — [Injul - C^(|ul)] - u Sj^(u) - cos u 
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g 2 (u) = u [Inlul - 


Cj(lul)l + f 


S^(u) + sin u 


and noting that Si(u) is an even function and gj(u) is an odd function, we get the 
distribution we shall call the “restricted energy loss distribution." 


exp[»c/l + 


f(X,A) = 


J cos[Xu + Kg 2 (u)Idu . (17) 


Comparing these results with the results of Vavilov's original derivation where 


f(X,A) 


f cos[Xu + Kg 2 (u)]|du 


X = - k(^» + 1 - t) 

^max 


8i(u) = /3*[ln|u| - Cj(|u|)l - u Sj(u) - cos u 


gjCu) = u[ln|u| - C^(Iu|) + ^*S^(u)l + sinu 


we note that has been replaced by f| in kX before the integral, and that X has an 

additional logarithmic term. Also. /?* has been replaced by eiP*/«max ^ exponential 
and in gi(u) and g] (u) . Thus, simply replacing «n,ax Vavilov’s distribution will 

not give the correct distribution. 
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To evaluate the sine and cosine integrals which occur in equations (1 S) and (16), we 
use the following approximations: 


t + In u 




G(u) sin u - H(u) cos u 


u £ 0. 01 

0. 01 < u < 1. 0 

1. 0 s u < 50. 0 
u SO.O 


® / 2j 1 

E (-l)^u •* 

(2J + 1)1 (2j * i) 

r/2 - G(u) cos u - H(u) sin u 

cos u 
ir/2 - 


u £ 0. 01 

0. OK u < 1. 0 

1. 0 £ u < 50. 0 
u 2 50. 0 


i a^u^ + a^ \ 

U \u“ + 850^ + EgU* + a^U'^ + Eg/ 

1 /u* + biU* + b,u* + b^u^ + bi \ 

tpVu'* + + bgu* + bju'* + 


3"' ; f,'^J 



» 38.027264 

as 

= 40.021433 

\ 

»2 

= 265.187033 

a« 

= 322.624911 


*» 

* 335.677320 

a? 

= 570.236280 

1 • ■ 

>4 

= 38.102495 

a. 

= 157.105423 






42.242855 


48.196927 


b. 

b, 

bj 

b4 

The series 


= 42.242855 
= 302.757865 
= 352.018498 
= 21.821899 


bs = 
b« = 482.485984 
bj = 1114.978885 
b, = 449.690326 


was terminated at the term which was 10'* of the sum to that term. 


Since most of the contribution in the integral in equation (17) is near the origin, the 
integration can be terminated at some finite limit. Simpson’s rule was used to perform the 
numerical integration. Difficulties are encountered in' obtaining good convergence of the 
integral where k is small and on the tails of the distribution where Ik | is large. Typically, 
for a (c of 0.87, good convergence was obtained on the tails [f(X,A) 10 *1 by 

integrating to a u of 20.0 with 6000 integration steps. At a k of 19.3, integration to 1.5 
with 2500 steps yielded good convergence. Good convergence here is taken as more than 
two significant figures since the test is applied at the extremes of distribution. Nearer the 
maximum of the distribution, convergence is much better than this. The test for good 
convergence was a second calculation where the upper bound of the integration was doubled 
and the number of integration steps was quadrupled. 


B. Results for the Restricted Energy-Loss Distribution 

In this section, we show some of the results obtained using the restricted energy-loss 
distribution derived above to calculate the energy-loss distribution due to low-energy 
transfers. 

Taking the ionization chamber in Reference 1 as a typical example, we show (Fig. 1) 
the energy-loss distribution for an iron nucleus incident at 2.5 GeV/nucleon on a 
5.1-cm-thick chamber containing 90-percent argon and 10-percent methane at 
1 atmosphere. The curve for “no cutoff” is the regular Vavilov distribution which as usual is 
skewed to the right, peaking at an energy slightly below the mean energy loss S of 
9.16 MeV. Tire stopping power S(E) used to calculate S was the proton stopping power 
at the same energy per nucleon multiplied by Z* . The electron density weight average of 
the argon and the methane stopping power tabulated in Reference 14 were used to 
approximate the stopping power of the mixture. As would be expected, as the cutoff e, is 
decreased, the tail to the right of the distribution maximum disappears, making the 
distribution more Gaussian and the mean of the distribution decrease. 


III. HIGH-ENERGY DELTA-RAY DISTRIBU'i'ION 


In this section, we describe the number and energy distribution of the high-energy 
delta rays to be used and test its validity. 
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Figure 1 . P.cstricted energy loss distribution with various delta-ray energy 
cutoffs for an iron nucleus (Fe*‘ ) with an energy of 2.5 GeV/nucleon for a 
pathlength of 5.1 cm in a mixture of 90-percent argon and 10-percent methane 
at 1 atmosphere of pressure, e—av was 12.853 MeV. 






Because equation (11) applies above the cutoff energy Ci , the average number of 
delta rays produced along the pathlength X above 6| is given by 
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( 20 ) 


Thus, the energy distribution of the high-energy delta rays is given by 




or 


n(e) = 


1 





€ € e 

max max max 


( 21 ) 


Since the production of a high-energy delta ray along the path of the primary 
particle path is a rather rare event, and since events of this type are independent of each 
other, we assume that the distribution of the number K of delta rays produced 
above e, is Poisson, with a mean of n^ or 


P(K) 



k: 


( 22 ) 


IS 
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To test the assumption of equation (22), we performed a small Monte-Carlo 
calculation to determine the distribution of energy lost to high-energy delta rays by primary 
particles. For each of a set of primary particles, we sampled a number K from 
distribution P(K) then we sampled K energies e from n(e) and summed them to obtain 
the energy loss of each primary particle. (See Appendix B for a description of sampling 
techniques.) Figure 2 shows the resulting distribution for 5000 iron nuclei incident on the 
argon-methane mixture, as in Figure 1 . The cutoff energy ei was 0.08 MeV. 



Figure 2. Energy loss distribution due to delta rays with energies above 0.80 MeV for an 
iron nucleus (Fe**) with an energy of 2.5 GeV/nucleon for a pathlength of 5. 1 cm in a 
mixture of 90-perccnt argon and 10-percent methane at I atmosphere of pressure. 


The total energy loss is the sum of the energy lost to low-energy delta rays and the 
energy lost to high-energy delta rays. According to Reference 1 2, if a ranciom variable Z is 
given by 

Z « t -t- u , 
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then the frequency distribution of Z is given by the convolution integral 


00 

fl(Z) = / ^ 2(2 - t) f3(t) dt 

.CO 


( 23 ) 


where fzCu) is the frequency distribution of u and fjft) is the frequency distribution 
of t . Since the distributions used here are zero over the negative range, the distribution of 
total energy loss f j(A j) is given by 


£t(At) = . t) dt (24) 


where Aj is the total energy loss and distribution of toss caused by 

high-energy delta rays, determined by the Monte-Carlo calculation. Figure 3 shows the result 
of applying equation (24) to the distribution shown in Figure 2 and the corresponding 
distribution in Figure 1. The distribution shown in Figure 3 is practically identical to the 
distribution in Figure 1 where no cutoff is used except in the range between 12.0 and 
13.0 MeV. (The discrepancies here are believed to be purely statistical in nature. The 
number of histories which gave contributions in this region was small.) Similar results have 
been obtained for other primary particle charges and energies and" for different energy 
cutoffs £i . Thus, we believe the assumption of equation (22) is valid. 


IV. HIGH-ENERoY. DELTA-RAY TRANSPORT AND ENERGY DEPOSITION 


In this section, we describe the method of transporting the high-energy delta rays 
produced in the ionization chamber and surrounding material and of calculating their energy 
deposition in the chamber gas. 

Since the chamber consists of a number of materials and has a fairly complex 
geometry, Monte-Carlo transport techniques are the only ones that are easily applicable. A 
pure simulation, following the consequences of each individual scatter as is sometimes done 
with neutrons or gamma rays, is not used, because an electron interacts with so many 
atomic electrons in the process of slowing down. Instead, the electron path is divided into 
several segments, and a new direction of travel is selected from a multiple scattering 
distribution at the end of each segment. The scheme used here follows closely those 
described by Berger [161 and Perkins [ 17] . First, an energy grid is set up uniformly spaced 
in the logarithm of energy between ei and ^niax’ energy cj is larger than the 

proceeding energy e|.| by a fixed percentage (usually about 5 percent). Second, using the 
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Figure 3. Total energy loss distribution calculated using the convolution j 

integral [equation (24)] on the distribution shown in Figure 2 and the ■ 

corresponding distribution (e I = 0.08 MeV) from Figure 1. ; 
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continuous-slowing-down approximation, the distance traveled AXj in losing 
energy e| - cj.i is determined from the total electron stopping power Syfe) [18] by 



/ 


de 

S^(^) 


(25) 


where p is the material density. Since £| is near £j.|, equation (25) is approximated by 


In 


AXj 


I-l 




w 


(26) 


which assumes a power law connecting points S-j<cj) and Sj(ej.]) on the stopping power 

curve. Where mixtures are involved, the stopping power is the average of the stopping 
powers of the elements or compounds in the mixture weighted by the corresponding 
electron densities. Third, given the energy lost and the distance traveled, the Moliere 
multiple scattering parameters (19,20,21) are calculated. (Appendix C describes the 
Moliere distribution.) 


The Moliere distribution was selected over other distributions that are theoretically 
more accurate because it can be very rapidly sampled. (The comparisons in Figures 4, 5, and 
6 show that the Moliere distribution is adequate for our purposes.) Since each primary 
particle can produce several hundred high-energy delta rays in passing through the ionization 
chamber and surrounding material, and since several thousand primary particles should be 
sampled to obtain satisfactory statistical accuracy, a multiple scattering distribution that can 
be rapidly evaluated is of prime importance in making the calculation possible. 

An electron history proceeds as follows: The electron initial energy is sampled from 
equation (21) (see Appendix B) and the nearest grid energy cj is taken as the actual 

starting energy (Fig. 7). The electron’s initial spatial coordinates are sampled randomly 
along the segment of the primary particle’s path in a region of a single material; i.e., if 100 
electrons are generated along a path segment through the chamber wall, their initial 
coordinates are distributed randomly along that segment. If five electrons are generated on 
the chamber segment, their initial coordinates are random along that segment. Energy 
deposition by the electron is determined by taking the difference Between the two possible 
sources of energy, the production of the electron in the chamber gas or the scattering of the 
electron into the chamber gas, and the one possible energy-loss mechanism, the escape of 


( 
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Z/r, FRACTION OF ELECTRON PATHLENOTH 

Figure 4. Energy deposition for 0.5 MeV electron normally 
incident on an aluminum slab. 


the electron from the chamber gas. If the electron initial coordinates lie in the chamber gas, 
the initial energy is added to the sum of energy deposited by the electron in that chamber. 
Otherwise, the sum is set to zero. If the electron escapes from the chamber, the energy at 
escape will be subtracted from the sum. The angle between the initial electron direction and 
direction of the primary particle is calculated from 


cos 0 ' = 




P 


fZ 


2mc^ 


(27) 
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Figure 7. Geometry' for high energy electron production and scattering 
showing electrons produced inside and exterior to the chamber. 
(The angles and lengths are not to scale.) 


[see Appendix A for a derivation of equation (27)1 ; the direction cosines of the initial 
electron trajectory are calculated using the cosines of the primary particle path (see 
Appendix D). Then, the coordinates of the end of the first path segment arc determined. 
Tests are made to determine if the chamber boundary has been crossed. If the electron has 
escaped from the chamber gas, the energy at the boundary determined by linear 
interpolation is subtracted from the energy sum deposited by the electron in that chamber. 
If the electron has scattered into the chamber gas, then the energy at the boundary is added 
to that chamber sum. When any cnamber boundary is crossed, a new set of segment end 
coordinates must be calculated because of changes in density or stopping power. Then, tests 
are made for succeeding boundary crossings until the region is reached where tne segment 
actually ends. A new direction is then sampled from the Moliere multiple scattering 
distribution, and tne process of boundary testing starts over along tne new path segment. 
The electron is followed in this fashion until it either escapes from the instrument or 
reaches 6i . At 6| , one last segment is taken whose length is the mean electron 
pathlength (181 at ci • 

As a check on the accuracy of the electron transport, calcula'.ions of energy 
deposition in a plane aluminum slab were made using the described method. Reference 22 
describes similar calculations using Berger’s program, ETRAN4161, which is recognized as 
one of the best electron transport codes presently available. Figures 4, 5, and 6 show 
comparisons of the results from Berger’s program [22] and the transport described here for 
monoenergetic electrons normally incident on the slab surface at 0.5, 1.0, and 4.0 MeV, 
respectively. The energy grid separation was 5.0 percent, and the cutoff energy ei was 
0.08 MeV in our case. Results from our program tend to peak slightly higher in energy and 
drop more rapidly with energy. The differences are small and should not significantly affect 
the results of our calculations. 


V. TOTAL ENERGY DEPOSITION 


In this section, we show how the total energy deposition distribution for primary 
particles is calculated using the results of the previous sections. 

In these calculations, the primary cosmic rays are presumed to arrive as a 
monoenergetic, monodirectional beam of a single isotope and to travel in a straight line 
through the instrument. Given the exterior point of incidence and the beam direction 
cosines, the pathlength in each chamber is determined and the low-energy, delta-ray energy 
deposition (restricted energy loss) distributions are calculated according to equation (17) for 
each chamber. 

To obtain the high-energy, delta-energy distribution, we use Monte-Carlo techniques 
in one of two schemes, depending on the number of delta rays produced by the primary. If 
each primary on the average generates only a few (< 100) delta rays above the cutoff 
energy, a new set of electron histories can be generated for each primary particle as the 
primary is followed through the instrument. (The number and energy sampling is as 
described in Section III.) In most cases of interest, however, each primary produces on the 
order of 1 000 high-energy electrons, and generating a new set of electrons for each primary 
particle becomes prohibitively expensive in terms of computer time. Instead, a set of 
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electrons is sampled from a table of electron histories previously generated and placed on 
magnetic tape, and the same electron history can be used several times. Tliis approximation 
is valid if the table can contain a representative population of electrons. For the calculations 
shown here, 20 000 electron histories were generated in each region of the geometry that 
was crossed by the primary beam. 

The electrons sampled for each primary particle are treated as follows: The initial 
energies of electrons that happen to be generated in a given chamber are summed and added 
to an energy loss/deposition for low-energy transfers sampled from the distribution given by 
equation (17) to obtain a total energy loss for the particle in the chamber. Next, the energy 
deposited in a given chamber by all electrons sampled for a primary is summed to obtain the 
energy deposited by high-energy electrons in the chamber. Combining this energy with the 
sampled low-energy electron loss/deposition sample yields the total energy deposited by the 
chamber primary. After several thous.and primary particles have been sampled, the resulting 
energy losses and energy depositions are sorted into energy bins to give histograms 
approximating the frequency distributions for low-energy ioss/deposition, high-energy 
deposition, total energy deposition, and total energy loss. 

A second and better method for calculating the total deposition 
distribution - because it eliminates the statistical fluctuation of sampling from the 
low-energy loss/deposition distribution — is to use equation (23) in a convolution of the 
restricted energy-loss distribution from equation (17) and the histogram frequency 
distribution for high-energy electron energy deposition resulting from the Monte-Carlo 
calculation. This is the method that is used for the results that follow. 


VI. RESULTS AND CONCLUSIONS 


Figures 8 and 9 show the total energy deposition distribution calculated by the pure 
Monte-Carlo (the histogram) and the Monte-Carlo convolution (the smooth curve) for iron 
nuclei Fe’‘ normally incident at 2.5 GeV/nucleon on the double ionization chamber of 
Israel et al. [ 1 1 . The instrument was simulated by a stack of five rectangular slabs 99-cm 
tong by 99-cm wide. The front slab which represents material between the primary beam 
and the first chamber plus the front chamber wall is 0.40 cm or 1.1 g/cm^ of aluminum 
thick. The second slab which reprt:ents the first chamber gas is 5. 1 -cm thick. The gas is a 
mixture of 90-percent argon plus 10-percent methane at I atmosphere with a density of 
1.6767 X 10'* g/cm^. Thus, 8.5 X 10'* g/cm’ of chamber gas is encountered by the 
beam. (The gas slab is identical to the one used for Figure 1 .) The third slab representing the 
back wall of the front chamber plus the front wall of the back chamber is 6.0 X 10"* cm or 
1.6 X 10*‘ g/cm’ of aluminum thick. The fourth slab is the gas of the second chamber and 
is identical to the second slab. The fifth slab is the back wall of the back chamber which is 
3.0 X 10’’ cm or 8.0 X 10 * g/cm* of aluminum thick. We assumed that little energy is 
deposited in the chamber by electrons reflected from the side walls of the chamber or from 
material behind the back chamber wall so that electrons that entered these regions are 
allowed to escape from the geometry. The cutoff energy e, was chosen to be 0.08 MeV so 
that comparisons can be made with curves in Figures 1, 2, and 3. 
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ENERGY DEPOSITION (M«Vl 

Figure 8. Calculated frequency distribution of energy deposition in the first of two stacked 
ionization counters of Israel et al. [ 1 ] by iron nuclei (Fe®‘) normally incident at 
2.5 GeV/nucleon. (The histogram is the pure Monte-Carlo calculation and 
the curve is the Monte-Carlo convolution.) 


Tlie distributions of Figures 8 and 9 appear to be Gaussian with no Landau tail 
present. The full width at half maximum is 0.9 MeV or about 10 percent of the mean energy 
deposition. The mean energy deposition is below the mean energy loss by about 3 percent, 
implying that energy deposition is not in equilibrium with energy loss. More energy is being 
transported out the back wall of the chamber than is being transported in the front wall, in 
agreement with Reference 23 which suggests that perhaps 6 g/cm* or more of material is 
required to reach equilibrium between the production and absorption of knock-on 
electrons. The most striking observation is the reduction of the full width at half maximum 
by a factor of two as compared with the width of the energy-loss distribution (Fig. 1 ). 
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Figure 9. Calculated frequency distribution of energy deposition in the second of two 
stacked ionization counters of Israel, et al. [ 1 J by iron nuclei (Fe'‘) normally 
incident at 2.5 GeV/nucleon. The histogram is the pure Monte-Carlo 
calculation and the curve is the Monte-Carlo convolution.) 



Observations by Israel et al. tl 1 on balloon flights from Palestine, Texas (where the 
geomagnetic cutoff is 4.5 MV), using the instrument simulated for Figures 8 and 9, indicate 
a width of 14 to 16 percent for the iron peak. (Their results are shown as the histograms in 
Figure 1 0 taken from Reference 1 .) The measurements had several sources of fluctuations 
not considered here. Amplifier noise was Gaussian with aj^ equal to 2.2 percent of the 

mean pulse height due to normally incident relativistic iron. Uricertainties in the trajectory 
measurement were characterized by a of 2.1 percent. The experimenters also suggested 
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Figure 10. Comparison of Israel et al. H ] results with our calculated 
distribution including noise and contribution from charges 25 and 27 







that part of the width of the peak was because perhaps as many as 10 percent of the counts 
included under the peak were for charges of 25 and 27. The superimposed curves in 
Figure 10 are out-results with the other sources of Huctuations folded into our calculated 
distribution and with a contribution of 10 percent equally divided between charges 25 and 
27. The curve is normalized to the total number of counts in channels 1 1 0 to 1 30, and 
energy deposition is assumed to be proportional to charge squared. The peak width with the 
other sources of fluctuation included is 1 3 percent of the mean energy deposition which is 
in reasonable agreement with their results. 

Figure 1 1 demonstrates the importance of the surrounding material for the energy 
deposition calculation and shows the energy deposition distribution for the bare 5.1 gas slab 
used for the previous figures with no surrounding material. The peak energy has de''reased 
by more than 1 .0 MeV and the full width at half maximum is nearly halved compared with 
the complete geometry. Note that the distribution is merely shifted slightly from the 
corresponding energy loss/deposition distribution for C| equal to 0.08 MeV in Figure 1 
implying that high-energy electrons from material above the chamber contribute more to 
energy deposition in the chamber than do those high-energy electrons created in the 
chamber. 

Figures 1 2 through 23 show the results of a study of the response as a function of 
charge of a proposed ionization chamber to be flown on a satellite [24). The instrument 
actually consists of four chambers with the responses of the first and second chambers 
combined to obtain a single measurement of energy depositionM The third and fourth 
chambers were similarly coupled. In the calculations, the instrument was simulated by a 
stack of 50- by 50-cm rectangular slabs. The first slab representing material in front of the 
instrument was 0.37 cm or 1 .0 g/cm* of aluminum. Next, theie were eight alternating slabs 
0.0508-cm or 0.4-g/cm* thick of iron representing the chamber walls and 3.0-cm thick 
representing gas. Last, there was the back wall of the last chamber which was like the other 
iron slabs. The gas was a mixture of 98-percent xenon and 2-percent carbon dioxide at 
1 atmosphere. Thus, the gas slabs were 1.73 X 10 * g/cm* thick. The cosmic rays were 
taken as normally incident at the center of the front aluminum slab with an energy of 
10.0 GeV/nucleon. 

For lower charges (Figs. 12 through 15), we note that the Landau tail is still in 
evidence but that the distributions are not as skewed as the Vavilov energy-loss distribution 
would be. The variation of the distribution maximum with charge can be well represented as 
being proportional to charge squared, and the variation of the distribution width seems to 
be proportional to charge. The width again is much smaller for the energy deposition 
distribution than it would be for the energy-loss distribution. 

In conclusion, the relationship between energy loss and energy deposited by primary 
cosmic rays is not simple. To determine the energy deposited in a small region, as is required 
for the analysis of data on heavy cosmic rays from ionization chambers and in other 
applications such as microdosimetry, the energy transport and deposition by high-energy 
secondary electrons must be considered. In the case of ionization chambers, the calculations 
described here show that the energy deposition frequency distribution is considerably 
nanower than the energy-loss distribution with the most probable energy deposition 
somewhat below the mean energy loss for the geometries considered. The deposition 
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Figure 1 1. Tlie energy deposition frequency distribution for iron (Fe*‘) 
nuclei normally incident at 2.5 GeV/nucieon on a bare slab with 
dimensions 99 X 99 X 5.1 cm composed of 90-percent argon and 
10-percent methane at a pressure of 1 atmosphere. 













I Figure 13. The energy deposition frequency distribution for 

I nitrogen (N'*) nuclei normally incident at 10.0 GeV/nucleon 

I on the four-chamber instrument. (These results are for 

I the second two-chamber-coupled detector.) 


32 











FREQUENCY 



ENERGY DEPOSITION IM«V) 


Figure 1 6. The energy deposition frequency distribution for sulfur 
(S’*) nuclei normally incident at 10.0 GeV/nucleon or the 
four-chamber instrument. (These results are for 

the first two-chamber-coupled detector.) -4 
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Figure 1 7. The energy deposition frequency distribution for sulfur 
(S^*) nuclei normally incident at 10.0 GeV/nucleon on the 
four-chamber instrument. (These results are for 
the second two-chamber-coupled detector.) 
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Figure 1 8. The energy deposition frequency distribution for chromium 
(Cr** ) nuclei normally incident at 1 0.0 GeV/nucleon on 
the four-chamber instrument. (Tliese results are for the 
first two-chamber-coupled detector.) 
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Figure 19. The energy deposition frequency distribution for chromium 
(Cr®*) nuclei normally incident at 10.0 GeV/nucleon on the 
four-chamber instrument. (These results are for the 
second two-chamber-coupled detector.) 
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Figure 20. The energy deposition frequency distribution for manganese 
(Mn**) nuclei normally incident at 10.0 GcV/nucleon on the 
four-chamber instrument. (These results are for the 
first two-chamber-coupled detector.) 
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Figure 21. The energy deposition frequency distribution for manganese 
(Mn**) nudei normally inddent at 10.0 GeV/nucleon on the 
four-chamber instrument. (These results are for the 
second two-chamber-coupled detector.) 
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distribution is narrower because the high-energy delta rays basically tend to follow the 
primary particle's path through the instrument. First, they are initially produced at small 
angles from the direction of the primary path, and, second, the multiple scattering 
distribution for high-energy electrons is highly peaked in the forward direction. As they 
travel through the instrument, they deposit their initial energies over long segments of the 
primary particle’s path, whereas in their production they have produced large fluctuations in 
the energy loss along short segments of the path. Because high-energy electrons produced 
above the chamber gas deposit more energy in the chamber gas than those high-energy 
electroiis produced in the gas, reductions of the matter above the chamber should decrease 
tlie width of the energy deposition frequency distribution and slightly lower its mean energy 
deposition which could improve the instrument’s charge resolution. 

Presently, the calculation of energy deposition for charges around 20 and above are 
almost prohibitively expensive in terms of computer time because of the large number of 
delta rays involved. With careful reprogramming, the computer time can probably be 
reduced. For example, electrons cannot travel farther than their mean pathlength, assuming 
no straggling. Thus, electrons produced one pathlength from the chamber need not be 
considered in the calculation. The second improvement that needs to be introduced for 
extension of the calculation to very high charges is the introduction of the correct delta-ray 
production cross section from a partial wave calculation [3] . With these improvements, the 
method could cover the whole charge spectrum. 
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APPENDIX A. DERIVATION OF THE ANGLE - ENERGY 
RELATIONSHIP FOR DELTA-RAY PRODUCTION AND 
THE EQUATION FOR MAXIMUM TRANSFERABLE ENERGY 


The initial relativistic momentum of the primary particle is given by 


P = yMV 


where V is its velocity, 7 is given by 


1-/32 


and M is its mass. Conservation of momentum gives 


P = P» + P 


where P' is the primary particle momentum and Pg is *he delta-ray momentum. Assuming 

that the primary particle was initially traveling along the X axis in the positive X direction, 
we may write equation (A-3) in component form as 


P = P* cos <t> +P cos 6 
e 


0 = P' sin (j> +P sin 0 
e 


where ^ is the polar angle of P’ and 0 is the polar angle of Pg . Since delta rays are less 
massive than nuclei, cos 0 is always greater than zero, and we may use the identity 
cos* 0 -► sin* ^ = 1 together with equation (A-5) to obtain 


cos 0 =*/l - sin 2 Q 

V p-« 
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Thus, equation (A-4) becomes 


P = P'ii I i ^ sin® 0 + P cos 0 


Using algebra and the identity cos* 6 + sin* d = 1 , we get 


P® - P'® + P® 


cos 6 = 


Using conservation of energy momentum 


E® = P? c® + M®. c* 
i 1 1 


and the momentum-energy relationship 


p 

i ~ c 


where Ej is the total energy, equation (A-8) becomes 


cos 0 = 


- E'® + E® - m® c* 
e 

2E/We® - m® c^’ 
e 


wJiere E is the initial primary particle total energy, E' is the final primary particle total 
energy. Eg is the delta-ray total energy, and m is the electron mass. Using conservation of 

energy 


E = E' + E - me® , 
e 
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equation (A-9) becomes 


cos 0 = 


E + me* 


E - me* 
e 


E + me' 
e 


(A-iO) 


Using the mass energy relationship 


E = yMc* 


and the definition of kinetic energy 


€ = E - me' , 
e 


equation (A- 10) becomes 


cos 0 = 


1 + M- 


fs V € + 2mc^ 


(A-11) 


For m/M ■< 1 , equation (A-1 1) is approximately 


e ^ ^ 777 ^ me* 


(A-1 2) 


Solving equation (A-1 1 ) for e we get 


2 me* fl* cos*0 




/3* cos* e 


(A-1 3) 


The maximum energy transfer occurs when 6 is zero; thus. 
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or squaring the bottom term and simplifying 



where m/M < 1 

€ fs 2mc® fl* y* 
max 






APPENDIX B. SAMPLING FROM FREQUENCY DISTRIBUTION 


A continuous frequency distribution f(X) , often called the density function, is 
defined such that f(X)dX gives the probability of an observation between X and 
X + dX . Tlie cumulative frequency distribution is defined as 


X 

F(X) = J f(X') dX' (B-1) 


It gives the probability of finding an observation X' less than or equal to X . For discrete 
frequency distributions, the cumulative distribution is given by 

N 

F(X^) = 2] P(XJ (B-2) 

1“0 * 


where P(Xj) is the probability of observing Xj . Obviously, F(X) varies between zero and 
unity. 

All sampling in these calculations is based on one simple rule 124): To obtain a 
sample whose frequency distribution approximates a given continuou' frequency 
distribution, the equation 

F(X^} = R (B-3) 


is solved for X|j the sampled variable where R is a pseudorandom number uniformly 
distributed on the interval 0 to 1 . For a discrete distribution, the inequality 




(B-1) 


is solved for Xj^fj . 

The pseudorandom number R is generated using the multiplicative congruence 
method. Given ine last pseudorandom interger Nj between 1 and 2*’ - 1 , the next 
integer in the sequence is given by 
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Nj + 1 = 5‘® Nj, mod 2” 


I + 1 
■^35 


As an example of sampling from a continuous distribution, consider tiie frequency 
distribution of the energy of high-energy delta rays given in equation (21). The cumula*:ve 
distribution is given by 


1_ -1_ + r In ^ 

„ , , ^max 

" ‘V- 1..1. 

Cl e e e 

* max max max 


Thus, we have from equation (B-3) 




€ C 

max max 


which can be solved numerically for by the simple iteration 


■[Jr 


+ — In — ^ 

€ Ct,! 

max Ri_l 


R ± -- i- . in -^1 

€i € e c 

* ma,x max max J 


Where the initial guess is taken as e, , convergence is rapid except at e^ax- 
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APPENDIX C. THE MOLIERE MULTIPLE SCATTERING DISTRIBUTION 


The Moliere muHiple scattering distribution gives the angular distribution expected 
for a beam of electrons which has been scattered several times (>20) . Bethe [21 J provides 
an easily understood description of the distribution. It is given by 


M(0} 0 do = (2e + g + + . . .) 


(C-1) 


where M(0) is the frequency of scattering through the angle 0 . The series is sufficiently 
accurate when terminated with the terms shown [21]. The 
functions M^’^(0) and are tabulated by Bethe [21]. The parameter B is 

defined by the transcendental equation 

B - In B = b (C-2) 


wl £re b is given by 


b = In 


1.16Tx„ 


Xg is the screening angle given by 


0.885a Z P 

o o e 


1.13 +3.76 



o 

137 iS 

e 


(C-3) 


where /i is Planck’s constant divided by 2 jt, a^ is the Bohr radius, Zq is the atomic number 
of the material, is the magnitude of the electron’s momentum, and (3g is the magnitude 
of the electron’s speed over the speed of light. Xg is given by 


_ 2e^ 

“ P /3 c 
e e 


^ ttN AX,Z (Z + 1}‘ 
a 1 o o 


(C-4) 
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"3 

! 

i 


where N„ is that atomic density (atoms/cm^ ). ’I' is defined by 
<1 


Since the electron energy changes in passing through AXj , the mean energy 
(cj + ej.j)/2 is used to define Pg and . Equation (C-2) is solved numerically for B by 

the simple iteration 


B. . = b + In B. 
i + 1 1 


The initial guess Bi for the initial scatter is 8.0. For succeeding scatters, the converged 
result from the previous scatter is used for the initial guess. Convergence to several places is 
usually obtained in five or six iterations of equation (C-6). 


Sampling from equation (C-1), we use equation (B-3) or 


f M(9) edO = R 


which reduces to 


A A 

r R r R 


/ ( 1 ^ 

R o M («) 
e + + 


/ 

o ('i') 

B^ 


= R (C-8) 


where, from equation (C-5), 
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The integrals in equation (C-8) are performed numerically to the tabulated values given in 
Reference 21. Since the exponential is the most important term on the right of equation 
(C-8), an initial guess at 'I' is obtained by solving the equation 




Then, the correct value is found by searching through the tabulated values of the integral 
from . Linear interpolation is used between tabulated values. 
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APPENDIX D. USE OF DIRECTION COSINES 
FOR DESCRIBING PARTICLE SCATTERING 


One of the most efficient ways of describing particle’s direction for the Monte-Carlo 
transport is the use of direction cosines. In Cartesian coordinates, a unit vector along the 
particle path is given by 

k" = 0(1 + yk (D-1) 


where a, p, and y are the X-, Y-, and Z-direction cosines, respectively. The utility of 
this method is shown by the calculation of the new particle coordinates when the particle 
has moved a distance r . The new coordinates are gi/en by 


X, = X + ar 
' o 

Yi = + ,3r 

Zt = + Vr 


where Xq, Y^, and Zq are the old coordinates. The problem arises in determining the 
new direction cosines when a particle scatters at some angle 6' from the original direction. 

If we define a new coordinate system with its positive Z axis directed along k' and 
its X axis directed along i’ where i' is defined by 


.. _ k' X k 

^ ~ Ik'xkl 


Since 


/i j 

k' X k = I a p 
\0 0 



(D-2) 


equation (D-2) becomes 
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o* + /3* 


‘ - I 2 = ^ ^ 
V a* + ^ 


For a right-hand coordinate system. 


I- _ X 

^ ~ Ik- X i' 1 


Since k' and i' are perpendicular, the m^gmtude of the cross product is one and 

j- = -=====, layi + ^-yj - (ot* + /J*) kl 

V a* + 

Now if we have a unit vector in the new system given by 


n = a-i' + /3-j- + y-k- 


it is g.vcn in the unprimed system by 


n = a- 


la^ + /3^' 


(fit - aj) 


+ /5- 




[cryi + ^yj - (a^ + 0^) kj 


+ y- [ai -I- iSj + ykj 


or regrouping 
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n 


+ 


V'a I 


+ 


+ 


a*P -t- P* ay 
V 0£* +/3* ' 

■“'“tJlfT' * T'/J i 

va* + /j2 

~P'\I a* + /J^' + y'y k 


Thus, the direction cosines of n in the old system are 

. a'P+P'ay 
a* = ^===:± — + y a 

■Ja^ +p^ 

/3* = + y'p (D-6) 

y* = -p + y*y 

Wc note that the direction cosine in the new system a' , P' , y' are related to the scattered 
angles by 

a' = cos <p ' sin 6 ' 

p' = sin sin 0' (D-7) 

y* = cos 0' , 

where 0' is the azimuthal scattering angle and 0 ' is the polar scattering angle. 
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APPENDIX E. DERIVATION OF THE CONTINUITY EQUATION 
OF THE ENERGY-LOSS FREQUENCY DISTRIBUTION 


Let ^ (XA) be the flux of particles that having passed through a distance X of 
material have lost energy A and let P(AX, be the probability that a particle with 

energy loss A' at X will have an energy loss A after passing through an additional 
distance A X of material where only interactions that result in energy losses below e i are 
considered to contribute to total energy losses A or A' . Then the flux 0(X + AX,A) is 
given by 


A 

<>(X+ AX,A> = <#» (X.A) + f P(AX,A'— A) <#) (X,A')dA' 

0 


00 

- / P(AX,A— A") <^(X,A) dA" 
A 


(E-1) 


where the first integral on the right is because of particles scattered into A and the second 
integral is because of particles scattered out of A . Since 


a<f» (XA) ^ lim d> (X 4- AX, A) - d> (X,A) 
ax AX--0 AX 


using equation (E-1 ) we obtain 


dd> (X.A } 

ax 


lim 

AX— 0 



P(AX,A^— A) 
AX 


^ (X, A ' ) dA ' 



P (AX. A— A") 
AX 


(X,A> dA" 


(E-2) 


Since we are considering only interaction which results in losses below energy e | in the 
limit of small AX, P(AX,A'‘*A) is zero for A'->A>C| and we may change the limits on the 
integrals in equation (E-2) accordingly. Thus, 



I (X,A') dA' 


A + ej 

- / w(E,A" - A) <?>(X,A) dA" 
A 

where 

! A - €j A > ej 
0 A < €j 

Letting e = A - A' , e' = A" - A , and b = A - B, we obtain 
b 

^ ay*™ “ / w(E,€; (X,A -€>d€ 

8X 


- f w(E,c'j d> (X,A> dc' 
0 


(E-3) 


(E-4) 


The probability density function f(X,A) is given by 


f(X,A> = — . (E-5) 

/ <)(X^'>dA' 

0 


Thus, the continuity equation of f(X,A ) is found by dividing both sides of equation (E-4) 
E 

by J (X,A ')dA' , yielding 
0 



8f (X.A) 
dX 


b 

= j* w(E,£) f(X,A - f)d€ 

0 

€ 

- J w(E,€'> f(x,€'> d€' . (E-6) 

0 
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